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ABSTRACT 

We place a limit on the logarithmic slope of the luminous quasar luminosity function 
at z ~ 6 of (3 <; —3.0 (90%) using gravitational lensing constraints to build on the 
limit of j3 <; —3.3 (90%) derived from an analysis of the luminosity distribution (Fan 
et al. 2003). This tight constraint is obtained by noting that of the two quasars which 
are lensed by foreground galaxies, neither are multiply imaged. These observations 
are surprising if the luminosity function is steep because magnification bias results 
in an overabundance of multiply imaged relative to singly imaged lensed quasars. 
Our Bayesian analysis uses the a posteriori information regarding alignments with 
foreground galaxies of the two lensed quasars, and provides a constraint on (3 that 
is nearly independent of the uncertain evolution in the lens population. The results 
suggest that the bright end of the quasar luminosity function continues to flatten out 
to z ~ 6, as is observed between z ~ 3 and z ~ 5 (Fan et al. 2001). Provided that SDSS 
Jl 148-5251 at z = 6.37 is magnified by an intervening lens galaxy at z ~ 5 (White et 
al. 2003), we also show that the high lens redshift in this system implies a co-moving 
density of massive galaxies that is close to constant out to high redshift. This is in 
agreement with the lack of redshift evolution in the velocity function of dark-matter 
halos with velocity dispersions near 200km sec -1 as predicted by the Press-Schechter 
formalism. The combination of constraints on the quasar luminosity function and lens 
galaxy evolution are used to compute an improved estimate for the z ~ 6 multiple 
image lens fraction of ~ 1 — 3%. 

Key words: quasars: general - cosmology:gravitational lensing 



1 INTRODUCTION 

The quasar luminosity function (LF) is the most basic 
property of the quasar population. At low redshifts sev- 
eral decades of study have yielded a well defined optical 
quasar LF with powerlaw slopes at both the faint and bright 
ends that do not evolve out to z ~ 3 (e.g. Boyle, Shanks & 
Peterson 1988; Hartwick & Schade 1990; Pei 1995a; Boyle 
et al. 2000). At higher redshifts only very bright quasars 
are currently observable at the magnitude limit of large 
surveys, and recent evidence suggests that the slope of 
their LF is significantly shallower than observed at z 3 
(Schmidt, Schneider, Gunn 1995; Fan et al. 2001a). There 
are now six quasars known with redshifts z J> 5.8 (Fan 
et al. 2001b, 2003). These very high redshift quasars pro- 
vide important constraints for studies of structure formation 
(Turner 1990; Haiman & Loeb 2001) and reionisation (e.g. 
Madau, Haardt & Rees 1999; Wyithe & Loeb 2003). Deter- 
mination of their LF (Fan et al. 2003) is critical if we are to 
address these issues. 



As a population of sources for gravitational lensing 
the z ~ 6 quasars are unique, being the only sample 
where the gravitational lensing probability may be of or- 
der unity (Wyithe & Loeb 2002a, b; Comerford, Haiman & 
Schaye 2002). The high expected lensing rate arises through 
a large magnification bias which increases the fraction of 
gravitational fenses at a given flux level by drawing sources 
from the fainter, more numerous population into a flux lim- 
ited sample. As a result the fraction of quasars in a sample 
that are multiply imaged by gravitational lenses is sensitive 
to the slope of the LF. Conversely, the observed multiple 
image lens fraction may be used to limit the unknown slope 
of a LF. This exercise was undertaken by Fan et al. (2003). 
They presented likelihood functions for j3 given the absence 
of multiply imaged quasars in the z ~ 6 sample, and found 
the lack of lenses to be surprising at the ~ 90% level if 
ft ^> —3.5. This lensing constraint is consistent with their 
findings for f3 through direct analysis of the luminosity dis- 
tribution in the sample which yielded f3 <; —3.35 (90%). 
However there is a complication. The lens fraction is lin- 
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early related (nearly) to the efficiency of the lens population, 
which is proportional to the expectation value of the veloc- 
ity dispersion to the fourth power and to the space density 
of galaxies. Moreover, the lensing rate requires extrapola- 
tion of the local galaxy population to higher redshifts, and 
is also sensitive to cosmology (e.g. Kochanek 1996). Thus 
the unbiased lensing cross-section is quite uncertain. 

While high resolution imaging data (Fan 
et al. 2001b, 2003) shows none of the six z ~ 6 quasars to be 
multiply imaged, galaxies have been detected near the line of 
sight to two quasars (Shioya et al. 2002; White et al. 2003). 
This is puzzling because magnification bias should result in 
highly magnified multiply imaged sources being over rep- 
resented among a population of quasars whose images are 
located near foreground galaxies. The effect becomes larger 
as the quasar LF becomes steeper. In this paper we present 
a Bayesian analysis that employs information on a posteriori 
alignments of the two lensed quasars. This statistic is much 
less sensitive to the uncertainties in the lens cross-section 
than the fraction of multiply imaged quasars, and we show 
that it produces a tighter limit on f3. 

The paper is set out as follows. In § 2 we compute the 
gravitational lens cross-section in light of the recently mea- 
sured velocity function of galaxies, the probability of multi- 
ple imaging for different LFs and the limits on the quasar 
LF that result from the observed lack of multiply imaged 
quasars. We then discuss a Bayesian approach to comput- 
ing the lens fraction in § 3. In § 4 we describe the two high 
redshift quasars thought to be magnified by gravitational 
lensing. The probability of getting a multiply imaged quasar 
within a sub-sample of quasars observed to be near a lens 
galaxy is discussed in § 5. These probabilities are used to 
compute likelihood functions for the fraction of lensed singly 
imaged quasars, and to derive limits on /3. In § 6 we use the 
redshifts of the lens galaxies to constrain simple parametric 
models for the evolution of the lens galaxy population. We 
then combine these results with the limits on j3 to estimate 
the expected multiple image fraction for z ~ 6 quasars in 
§ 7. Finally in § 8 we discuss the implications of possible 
multiple imaging in SDSS J1148-5251 for /3 before present- 
ing our conclusions in § 9. Where required we assume the 
most recent cosmological parameters obtained through fits 
to WMAP data (Spergel et al. 2003). These include density 
parameters of Q. m = 0.27 in matter, Q b = 0.044 in baryons, 
SIa = 0.73 in a cosmological constant, and a Hubble con- 
stant of H = 71 kms -1 Mpc -1 . 



2 THE GRAVITATIONAL LENS 

CROSS-SECTION AND MULTIPLE 
IMAGING RATE 

The probability that a Spherical Singular Isothermal (SIS) 
galaxy will lens a back-ground source is proportional to the 
4th power of its velocity dispersion a. Thus what is re- 
quired to compute the cross-section for gravitational lens- 
ing is the velocity function. Until recently the velocity func- 
tion for early type galaxies (the dominant lens population, 
Kochanek 1996) had to be computed through combination 
of a galaxy LF with the Faber- Jackson (1976) relation. How- 
ever this procedure ignores the intrinsic scatter in the Faber- 
Jackson (1976) relation and is an unreliable method (e.g. 



Kochanek 1993). A more reliable representation is now pos- 
sible using the measured velocity dispersion function of early 
type galaxies. Sheth et al. (2003) presented the measured 
velocity dispersion function for early type galaxies. They 
suggested an analytic fit of the form 



(a)da - 



exp [— ((j/cr*)' 3 ] da 



r(a//3) 



P- 



(1) 



where <f> t is the number density of galaxies and a* is a 
characteristic velocity dispersion. Sheth et al. (2003) found 
that the parameters a*, a and (5 are strongly correlated 
with one another. From Sheth we take 0* = (2.0 ± 0.1) x 
10~ 3 (h^Mpc)- 3 , a = 6.5 ± 1.0, f3 = (14.75/a) ' 8 and 
a„ = 16ir(a//?)/r[(a + l)/f3] km s _1 . 

Given the Einstein Radius (ER) for an SIS 



Co = 4tv [ - 



(2) 



and a constant co-moving density of galaxies, the multiply 
imaged gravitational lens cross-section for sources at z s is 



dz 



da(l + z) 



cdt 2 

)-7-7r£(>- 

dz 



(3) 



We found the mean and twice the variance for a set of r com- 
puted assuming the parameters 0* and a to be distributed 
as Gaussian within their quoted uncertainties. This proce- 
dure yields r(z s = 6) = (2.5 ± 0.25) x 10~ 3 . The statistical 
uncertainty of 10% is significantly smaller than is obtained 
through use of the Faber- Jackson (1976) relation and a LF. 
The value of t(z b — 6) obtained from the velocity function is 
a factor of ~ 3 smaller than obtained in estimates of the lens 
fraction in the very high redshift quasar samples (Wyithe & 
Loeb 2002a,b; Comerford, Haiman & Schaye 2002). While 
the implied lens fraction for these samples is still expected to 
be an order of magnitude higher than in lower redshift sam- 
ples, it may result in less than one lens being expected in the 
current sample, hence limiting the use of the lens fraction 
for constraining the slope of the LF (Wyithe & Loeb 2002b; 
Comerford, Haiman & Schaye 2002; Fan et al. 2003). 

As shown above the use of a measured velocity func- 
tion reduces the statistical uncertainty in r to 10%. How- 
ever the remaining systematic error due to the uncertain 
evolution in the lens population (e.g. Keeton 2002) makes 
limits based on the lens fraction uncertain. For example, 
if the co-moving density of galaxies drops in proportion 
to (1 + z)~ 7 (density evolution) then for 7 = 1 we find 
t(z s = 6) = (0.9 ± 0.09) x 10~ 3 , while for 7 = 2 the value 
drops to t(z b = 6) = (0.4 ± 0.04) x 10~ 3 . Unless otherwise 
specified we assume a constant co-moving density of galax- 
ies (7 = 0). We consider only spherical lenses in this paper. 
Previous studies (Kochanek & Blandford 1987; Blandford & 
Kochanek 1987) have found that the introduction of ellip- 
ticities 5S 0.2 into nearly singular profiles has little effect on 
the lensing cross-section and image magnification. However 
the strong magnification bias will favour a high fraction of 
4- image lenses (Rusin & Tegmark 2001). Finally we note 
that the i-band dropout quasars are selected independent of 
morphology, and so do not select against lenses, though the 
possibility that the lens galaxy itself will prevent detection 
should be accounted for (Wyithe & Loeb 2002b). 

With these points in mind we compute the a-priori 
probability of multiple imaging in a sample of six quasars 
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Figure 1. Left: Likelihood functions for fi. Centre: Differential probability distributions for fi. Right: Cumulative probability distributions 
for fi. The lensing constraints based on the fraction of multiply imaged quasars in the sample are shown by the dashed grey curves. The 
lensing constraints that include information on the alignments between quasars and foreground galaxies are given for different values of 
x (dark lines). The solid grey curves in each panel correspond to the likelihood and probability functions for fi based on the luminosity 
distribution alone (Fan ct al. 2003). Results arc shown for a double powerlaw LF with a = —1.75 and AM = 4. A constant co-moving 
density of lens galaxies was assumed. 



with z ~ 6 for different LFs. A similar calculation has al- 
ready been performed by Fan et al. (2003), though we add 
the improvement of a more accurate r computed from the 
velocity function, extend the calculation to account for con- 
straints from the luminosity distribution (Fan et al. 2003), 
and provide a posteriori limits on f3. We also modify r to 
include only those lens galaxies which would not have con- 
taminated the i— band dropout selection of the quasars as 
discussed in Wyithe & Loeb (2002b). The probability that 
a quasar with luminosity L and magnification bias B will be 
multiply imaged by a foreground galaxy is 



Plena (L) 



B(L)t 



(4) 



B(L)t+(1-t)' 

where pi cns is slightly overestimated because we have not 
included any magnification bias for single image quasars. 
The sum of magnifications of multiple images (//) formed by 
an SIS has a probability distribution of the form 

dPm 8 



for fi > 2, 



(5) 



resulting in a magnification bias B(L) for a SIS and a LF 
$(L) of 



B(L) = 



rf/i 8 $(L/n) 
li M 3 *(L) ' 



(6) 



The probability given a LF <£>, that in a sample of six quasars 
at z ~ 6, none will be lensed is 



Pnolcns = [J (1 - Plena (ii)) , 



(7) 



i=l,JV q 



where JV q = 6 and the Li are the luminosities of the JV q 
quasars. 



2.1 lens fractions for double powerlaw luminosity 
functions 

A successful fit to the low redshift (z < 2) quasar LF is the 
double powerlaw (e.g. Boyle et al. 2000) 



*(L) 



$ 



(8) 



where a and (3 (note we have defined a and (3 to be nega- 
tive) are the slopes of the faint and bright ends of the quasar 



LF respectively, and we have expressed the luminosity L 
in units of the characteristic break luminosity. Use of the 
double power-law form for <J?(L) implies that we must spec- 
ify two additional parameters before deriving limits on (3. 
Firstly, since high magnifications will draw quasars that are 
fainter than the break into the sample, the magnification 
bias will depend on a. Moreover use of the double powerlaw 
requires specification of the quasar luminosity with respect 
to the break. The z ~ 6 quasars range in luminosity from 
M1450 = -27.15 -> -27.90 (Fan et al. 2001b, 2003). In this 
paper we specify the luminosity of the quasars relative to the 
LF in terms of the difference between the magnitude of the 
faintest quasar and that of the break (AM = 2.51og 10 L). 
For example if the break were at M1450 = —26.15 then 
AM = 1. 

The resulting probability p n oiens that in the sample of 
six z ~ 6 quasars, none will be multiply imaged (given a — 
— 1.75, AM = 4) is plotted as a function of /3 in the left 
hand panels of figure 1 (thick dashed grey line). The lack 
of multiply imaged quasars is only surprising at the 50% 
level for values of (3 < —3.3. The probability p n oiens is also 
tabulated in Table 1 for various values of a, (3 and AM. The 
lack of multiply imaged quasars is only surprising at the 10% 
level if (3 <, -4 and AM £ 3 or if (3 <, -3.75 and AM £ 4. 
We note that these probabilities (and those in figure 9 of 
Fan et al. 2003) represent a likelihood function for (3 rather 
than direct to limits on (3. We now turn to computation of 
these limits. 

The posterior probability for f3 is 



dP 



where 



Np n 



dP D1 



(9) 



dfi 



is the prior probability for (3, and N is a nor- 
malising constant. We assume that the prior probability for 
the slope is flat between two bounds /3mm and /3 m ax, hence 

dP v , 



d(3 



/3n 



(10) 



The absence of multiply imaged quasars in the sample is 
surprising if j3 is small because the multiple image magnifi- 
cation bias is large in that case, but is less surprising as (3 is 
increased. Hence the likelihood function p no icns is increasing 
with (3 so that the fraction of multiply imaged quasars car- 
ries no information on its upper limit. As a result, the con- 
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Table 1. The probability given /3, that in a sample of six z ~ 6 quasars none will be lensed. The values are tabulated assuming a constant 
co-moving density of lens galaxies and various values of a and AM. 





a 






1.0 








1.75 








2.5 






AM 


1 


2 


3 


4 


1 


2 


3 


4 


1 


2 


3 


4 


p 


= -4.0 


0.77 


0.42 


0.11 


0.01 


0.70 


0.35 


0.08 


0.006 


0.47 


0.18 


0.03 


0.002 


p 


= -3.75 


0.82 


0.57 


0.28 


0.08 


0.77 


0.51 


0.24 


0.07 


0.56 


0.33 


0.13 


0.03 


p 


= -3.5 


0.87 


0.70 


0.50 


0.32 


0.83 


0.67 


0.48 


0.29 


0.66 


0.51 


0.37 


0.23 


p 


= -3.0 


0.92 


0.86 


0.78 


0.72 


0.89 


0.84 


0.78 


0.72 


0.78 


0.75 


0.72 


0.70 


p 


= -2.5 


0.96 


0.93 


0.91 


0.90 


0.94 


0.93 


0.91 


0.90 


0.87 


0.87 


0.87 


0.87 


p 


= -2.0 


0.97 


0.96 


0.96 


0.96 


0.96 


0.96 


0.96 


0.96 


0.92 


0.91 


0.90 


0.89 



fidence with which a small value of can be excluded given 
the lensing observations alone is sensitive to /3 max . This de- 
pendence implies that a small probability p given some 
that a sample of quasars would contain no lenses does not 
translate to a limit on at the 100(1 —p)% level. A second, 
independent constraint that limits the upper bound on 



and hence lowers the dependence on 



d/3 



is required. 



2.2 the addition of constraints from the 
luminosity distribution 

Fan et al. (2003) have derived the constraint ± A/3 where 
= —2.3 and A/3 = 0.8 from their analysis of the luminosity 
distribution of the z ~ 6 quasars. We take this bound to be 
the 1-cr level of a Gaussian likelihood function for 0. We 
have plotted this likelihood function for /3 (normalised to a 
maximum of 1), as well as the differential and cumulative 
probability distributions for in the in the left-hand, central 
and right-hand panels of figure 1 (solid grey lines) . 

We now have two constraints on 0, one from the fraction 
of multiply imaged quasars and one from the distribution of 
luminosities. Assuming these constraints to be independent 
(this assumption is discussed below), we find a joint likeli- 
hood function, and hence posterior probability distribution 
for 



dP 

d0 



-/Vpnolens exp 



A/3 , 



dP p I'll ) 



, (11) 



The likelihood function is now normalisable because the dis- 
tribution of luminosities constrains large values of (Fan et 
al. 2003), and is quite insensitive to the prior ^ | N as 
a result. We plot the posterior differential and cumulative 
probability distributions in the central and right hand pan- 
els of figure 1. These are shown by the dashed light lines, 
and should be compared to the posterior probability distri- 
butions based on the distribution of luminosities alone (solid 
light lines) . The lensing constraint disfavours smaller values 
of 0, resulting in a narrower probability distribution for 0. 
The most likely value for the slope is ~ —2.2. The addition 
of the constraint from the fraction of multiply imaged lensed 
quasars improves on the limits obtained by Fan et al. (2003) 
from the distribution of luminosities alone. For this choice 
of AM and a, the lack of multiply imaged quasars rules out 
—3.0 at the 90% level. We also construct posterior cu- 
mulative probabilities for various values of a and AM and 
plot the results in figure 2 (thick dashed light lines). The 



addition of lensing constraints significantly improve the LF 
limits provided that the quasars are not too close to the 
LF break. For a constant co-moving density of lens galax- 
ies (upper two rows) we find <; —3.2 — > —3.0 (90%) for 
AM ;> 2. 

It should be noted that the two constraints are not quite 
independent. In general gravitational lensing tends to flatten 
the slope of the quasar LF by drawing populous faint quasars 
into a bright quasar sample (e.g. Pei 1995b). However there 
are two reasons to think that this is not a problem within the 
very high redshift quasar sample. Firstly the objects that are 
lensed are not multiply imaged and so have magnifications 
smaller than ~ 2. Secondly, the average change of slope is 
A/3 ~ 0.2 even in the most optimistic lensing scenario (where 
~ -3.5) for the z ~ 6 quasars (Wyithe & Loeb 2002b). 



3 BAYESIAN APPROACH TO COMPUTING 
THE MULTIPLE IMAGE FRACTION 

We may also use a Bayesian approach to compute the multi- 
ple image fraction. This will provide us with a natural frame- 
work within which to add additional a posteriori informa- 
tion on alignments of quasars with foreground galaxies in 
§ 5. Consider sources with unlensed impact parameters (in 
units of the ER) y = x — 1 with associated magnifications 
/x. We write the likelihood per logarithm of x of observing a 
singly imaged lensed quasar (including magnification bias) 



Lgingle — x{x 1)- 



1 $(L/^ sing l e ) 



/^single 



(12) 



where the factor (x — 1) accounts for the additional solid 
angle available at large y, and /i s i ng ie = x/(x — 1). This 
likelihood may be compared to the corresponding average 
likelihood of observing a multiply imaged quasar 



1 3> (L fi) 
dyy — ZJ7Y~ 



1 r *E± 

2 J 2 V V 



8 _ B{L) 



*(L) 



(13) 
(14) 



The likelihood that a quasar will be singly imaged at x 
rather than multiply imaged is therefore 



Psingle(/3 I X) = 



-^single 



-^single ~\~ -^mult 



(15) 



We may also calculate the posterior probability that a 
quasar will not be lensed 
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pnolcns,baycs 



s: 



dxp aing i e (l3\x) 



dP p 
dx 



(16) 



where — gp 1 is the prior probability for the x. This quantity 
is one minus the lens fraction and may be approximated 
using the usual formula for the lens fraction rB or more 
accurate forms such as equation (4). 

For the sample of z ~ 6 quasars, the posterior proba- 
bility distribution for (5 is therefore 



dP 

dp 



— N exp 




■ (17) 



where the prior probability distribution for x can be com- 
puted from the derivative of the Poisson probability that a 
source lies within a circle of radius x — 1 around a randomly 
positioned galaxy 



dx 



2t(x- l)e~ 



r(x-l) 2 



(18) 



Equation (17) yields identical limits to those based on the 
multiple image fraction (dashed light lines in figures 1 and 2) 
as computed in the usual way from equation (4). Moreover, 
the magnification distribution for singly imaged sources is 
naturally normalised within the formalism, and hence the 
magnification bias of singly imaged sources is directly in- 
cluded in the calculation. 



4 TWO LENSED Z ~ 6 QUASARS 

While none of the six z ~ 6 quasars discovered by Fan et 
al. (2001b,2003) in the Sloan Digital Sky Survey data are 
multiply imaged, two have close alignment with a foreground 
galaxy, implying that they are moderately magnified. We 
therefore refer to these quasars as lensed, though neither is 
multiply imaged. 

• SDSS J1044-0125 at z = 5.74: Shioya et al. (2003) 
have reported a faint foreground galaxy with z ~ 1.5 — 2.5 
at a separation of 9 = 1.9". Shioya et al. (2003) estimate 
the velocity dispersion to be a ~ 140 — 280km/sec for 
this redshift range. A second image would be detectable if 
a <; 220km/sec. For a SIS galaxy this implies that the mag- 
nification of the image could be as high as fi — 2, and that 
the image is located at x ~ 2 — 10 ER (where x ~ 10ER 
corresponds to a a ~ 140km/sec SIS at z ~ 2.5). 

• SDSS J1148-5251 at z = 6.37: This is one of 
two quasars known to exhibit Gunn-Peterson absorption 
troughs. White, Becker, Fan & Strauss (2003) present spec- 
tra showing emission features in the Ly f3 trough which they 
interpret as being Ly a emission from a foreground galaxy 
at z ~ 4.9. This is a likely scenario if gravitational lensing is 
important and was predicted by Wyithe & Loeb (2002b). 
Indeed, the smaller than expected Stromgren sphere im- 
plies that this quasar has been magnified by the intervening 
galaxy (White et al. 2003). While the alignment is presum- 
ably high, we do not know the degree of alignment with the 
foreground galaxy, and hence we cannot know x, though is 
the quasar is not multiply imaged it must be larger than 



2, and is probably comparable to the case SDSS J1044- 
0125. White et al. (2003) do however provide an estimate 
of a ~ 250km/sec for the foreground galaxy from the veloc- 
ity structure seen in the C IV absorption system. For an SIS 
at z — 4.9 lensing a source at z = 6.37 with a separation 9 
we have 



250km/sec 



ER. 



(19) 



5 IMPROVED LIMITS ON f3 FROM CLOSE 
ALIGNMENTS WITH FOREGROUND 
GALAXIES 

The Bayesian approach (§ 3) to computing limits on /3 from 
the multiple image fraction allows us to include the a pos- 
teriori information on the alignments of the two lensed the 
z ~ 6 quasars. 

For illustration, we begin with a hypothetical sample of 

6 quasars with AM — 4, a = —1.75 and (lensed) impact 
parameter x. The relative likelihoods for different values of 
P given the lack of multiple images are [p s in g ic(/3 | x)] Nq as 
specified in equation (15) where N q — 6. These likelihoods 
are shown in the left hand panel of figure 1 for values of 
x ranging from 5 to 20 (dark lines). Smaller values of (5 
are strongly disfavoured, particularly if x is not too large. 
Next we find find the joint likelihood function and hence 
a posterior probability distribution for /3 given a common 
impact parameter x for six quasars that are not multiply 
imaged 



dP 

d/3 



A^bsingiol^exp 



1 (p-p 




(20) 



The resulting posterior differential and cumulative proba- 
bility distributions are shown in the central and right hand 
panels of figure 1 (dark lines). The most likely value is 
near p ~ — 2 and the cumulative distributions suggest that 
P —3.1 — > —2.7 at the 90% level where the systematic 
dependence is on x. Thus any additional information about 
close alignments produces constraints that may significantly 
tighten the lower limits on the slope of the z ~ 6 quasar LF, 
with improvements in the limit that are greater than 0.5 
units in P for cases where the alignment is high. 

The above example suggests strong dependence of the 
limits derived for P on the value of x. For the two lensed 
SDSS quasars discussed in § 4, there are observational limits 
on x in the form 

dP _ ^ rf-Pprior (21) 

dx x dx ' 

where L x is the likelihood for x given the observations of a 



and z for the lens galaxy, and 



is the prior probability 



for x. Given the relation x = x(a, z) the likelihood L x is 

j j (^X Ox G?Pprior ^Pprior /oo\ 

x =LcrL z -— — . 

da oz da dz 

For SDSS J1044-0215 Shioya ct al. (2003) find the major- 
ity of the dependence in the likelihoods for a and z to be 
systematic, while there is no information in this regard for 
SDSS J1044-0215 (White et al. 2003). We assume flat dis- 
tributions ^ with limits of 2 < x < 10 for the two lensed 
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Figure 2. Cumulative probability distributions for j3. The lensing constraints based on the fraction of multiply imaged quasars in the 
sample are shown by the dashed grey curves. The lensing constraints that include information on the alignments between quasars and 
foreground galaxies are denoted by the dark lines. The solid grey curves correspond to the probability functions for /3 based on the 
luminosity distribution alone (Fan et al. 2003). Results are shown for a double powerlaw LF for various values of a and AM. The upper 
and lower two rows show results assuming density evolution with 7 = and 7 = 2 respectively. 



quasars. For the other 4 quasars we assume the prior prob- 
ability distribution for x (equation 18). 

The posterior differential probability distribution for /3 
then becomes 



dP AT 



n 




/ ) d(3 

dxip Binglc (P\ x,Li) 



n 



dP rior 

efep s in g ie(/3| x, Li) *j£ 1 



(23) 



The resulting cumulative probability distributions are 
shown in figure 2. The limits on f3 are significantly tighter 
than those obtained from the distribution of luminosities 
alone except in cases where a is large (shallow faint end 
slope) and AM is small (so that magnifications associated 
with multiple images tend to draw quasars with AM < 
into the sample) . In addition the limits are tighter than those 
obtained through consideration of the luminosity distribu- 
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tion and multiple image fraction (dashed light lines). For a 
constant co-moving density of lens galaxies and AM <; 2 
we find (3 ^ -3.1 -> -3.0, while for AM ^ 4 we obtain 
/3 k, -3.0 -> -2.9 (both with 90% confidence). Thus the 
tightest limits come from the inclusion of the a posteriori 
information that two of the quasars have close alignment 
with foreground galaxies. 

A second important point regarding the limit on /? pro- 
vided by equation (23) is that unlike the limit from the mul- 
tiple image fraction, it is nearly independent of the value of 
r. To demonstrate this independence we have computed con- 
straints on /3 (shown in the lower two panels of figure 2) that 
assume a dependence in the co-moving density of galaxies of 
(l + z) -7 where 7 = 2 (resulting in r = 0.0004). These limits 
may be compared with results that assume a constant co- 
moving density of lens galaxies (in the upper two panels of 
figure 2). The limits obtained from the multiple image frac- 
tion are much weaker if 7 = 2. On the other hand the limits 
that use a posteriori observations of the quasar-lens galaxy 
alignment are quite insensitive to 7. The reason is that the 
role of t is replaced by ^ for the two quasars which provide 
the largest contribution to the likelihood change between 
large and small values of /3. 



5.1 a posteriori choice of statistic 

We have computed limits on the value of (3 using two lens- 
ing based constraints, and a posteriori chosen the better 



one. This practice becomes unfair if a large number of dif- 
ferent constraints are available where each produces a dif- 
ferent limit. In the situation described we have two different 
lensing constraints. However the second constraint utilises 
additional rather than different information. Thus we are 
justified in choosing it a posteriori. 



6 LIMITS ON GALAXY EVOLUTION FROM 
LENS GALAXY REDSHIFTS 

As noted by White et al. (2003) the candidate lens galaxy 
in the system SDSS J1148-5251 is found at an improbably 
high redshift, which could provide an argument against the 
lens hypothesis in this system. This is quantified in figure 3, 
where the upper curves show the differential cross-section 
for a source at 2 = 6.37 lensed by a constant co-moving 
density of galaxies (left panel), and the fraction of the total 
cross-section that is found at a redshift larger than z (right 
panel). We see that the prior probability of finding a lens 
at z J; 4.94 among two lensed quasars is only ~ 0.01. How- 
ever the probability of finding a lens at high redshift is an 
a-posteriori statistic, i.e. we have chosen one of a possible 
number of a-priori unlikely events after the observation has 
already been made. Moreover the selection of lenses within 
the sample is not uniform in redshift. In particular, since the 
galaxy in front of SDSS J1148-5251 was identified spectro- 
scopically via its Ly-a emission line, it would be more easily 
identified at high redshift. In addition, the cross-sections 
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Figure 4. The velocity dispersion function of dark matter halos 
at several redshifts. 



plotted in figure 3 refer to multiple imaging so that the low 
probability for a high redshift lens results from the small size 
of the Einstein Ring radius at high redshift. In contrast, the 
redshift distribution of galaxies (with velocity dispersions 
larger than a) that lie within some specified angular separa- 
tion along the line of sight to a back-ground quasar implies 
that there would be 1 chance in ~ 3 of finding such an 
aligned galaxy at z > 4.96. The role of magnification bias 
in making lensing a likely scenario probably results in an a- 
posteriori probability of finding a high redshift fore-ground 
galaxy that lies somewhere between these two extremes. 

The redshift distribution of gravitational lenses may be 
used to constrain evolution in the lens galaxy population 
(Kochanek 1992; Ofek, Rix & Maoz 2003). While the ab- 
solute probability of having observed a high redshift lens is 
difficult to quantify, we may more easily discuss the rela- 
tive likelihoods of observing a high redshift lens as a func- 
tion of lens galaxy population. For definiteness we consider 
two parameterisations for the evolution of the lens galaxy 
population: firstly evolution of the characteristic velocity 
(Ti,(z) = <r*(l + z)~ 7 , and secondly evolution of the char- 
acteristic density (j>* = <^*(1 + z) ' , which we term veloc- 
ity and density evolution respectively. The left hand panels 
of figure 3 demonstrate the effect on the lens cross-section 
of varying 7. Values of 7 that differ from (constant co- 
moving evolution) result in a lens population that is trun- 
cated at high redshift. This effect may also be seen in the 
central panels of figure 3 where we have plotted the fraction 
of cross-section at redshifts larger than z. Values of 7 > are 
disfavoured by the existence of a lensing galaxy at z = 4.96 
among a sample of only two lens galaxies. 

To quantify this statement we construct a likelihood 
function for 7 from the product of the normalised probabili- 
ties for the lens redshifts. The likelihood should include con- 
stants (si) to account for the relative detectabilities of the 
two lenses (the two galaxies were discovered separately via 
different techniques) , though the limits on galaxy evolution 
are independent of these since the constants are independent 
of the evolution. The likelihood function is 



where the =L are differential cross-sections evaluated at the 

az 

lens and source redshifts, the integrals over the distribu- 
tions ^-{x — l) 2 account for the relative alignments of the 
quasar and galaxy, and the Bi are the magnification biases. 
The relative likelihood is dependent only on the product 
of the differential cross-sections. Note that the likelihood 
(equation 24) is only applicable if the magnification bias 
has aided in selection of the quasar, so that the source may 
be considered lensed. This is a caveat to the constraints im- 
posed on 7 in this section. However the quasar does appear 
to be magnified, as evidenced by the smaller than expected 
Stromgren sphere (White et al. 2003). 

In the right hand panels of figure 3 we plot the posterior 
cumulative probability for 7 



P(<7) 



-f 



AT T dPp T i OI 



(25) 



assuming a flat prior probability distribution for 7 at val- 
ues greater than 0. By excluding the possibility of 7 < 
we are assuming that the lens galaxy population increases 
monotonically in time as expected in hierarchical merging 
scenarios. This choice also leads to more conservative limits 
on 7. We find 7 0.4 and 7 1.6 at the 90% level assum- 
ing velocity and density evolution respectively. The possible 
presence of a lens galaxy at such a high redshift therefore 
offers an opportunity to constrain the (mass selected) co- 
moving density of massive galaxies to be close to constant 
out to high redshifts. This result is consistent with the study 
of Ofek, Rix & Maoz (2003) who performed a detailed study 
on a large sample (15) of multiple image lenses at z ~ 1 — 2. 

A constant co- moving density of lens galaxies out to z ~ 
5 may not be surprising in light of the Press-Schechter (1976) 
prediction for the velocity function of dark-matter halos 
(number per cubic comoving Mpc per unit velocity). Taking 
the circular velocity « v ir to equal the virial velocity of an 
SIS dark matter halo with mass M (Barkana & Loeb 2001), 
and a velocity dispersion a — v v i T /\/2 we find 



dn 
da 



— — M 



dn \/2 

dM 3«vir 



(26) 



where is the Press-Schechter (1976) mass function. The 
resulting velocity function of dark matter halos is plotted 
in figure 4 at a series of redshifts. Note that near velocity 
dispersions of a ~ 200km sec - , which dominate the lens 
cross-section, there is little evolution in ^ from z ~ 1 to 
z ~ 5. If massive galaxies occupied dark matter halos in 
the past as they do today, we would therefore expect little 
evolution in the lens population, even out to large redshifts. 



7 WHAT IS THE MULTIPLE IMAGE 

LENSING RATE FOR Z ~ 6 QUASARS? 

We may combine the information obtained for (3 and 7 
and estimate the expected multiple imaging rate for z ~ 6 
quasars. Since our earlier limits on f3 are independent of r, 
while the limits derived for 7 are insensitive to the magnifica- 
tion bias, the constraints placed on (3 and 7 in § 5 and § 6 are 
independent. In figure 5 we have plotted the resulting joint 
probability function (dark contours); the contours shown are 
at 61, 26, 14 and 3.6% of the peak value, corresponding to 
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Figure 5. Dark contours: the joint probability function for (i and 7. The contours show 61, 26, 14 and 3.6% of the peak value, 
corresponding to the 1, 2, 3 and 4 a levels. Grey contours: The predicted lensing rate. The solid, dashed, dot-dashed and dotted 
contours correspond to lens fractions of 0.1, 0.03, 0.01 and 0.003. The upper and lower rows correspond to velocity and density evolution 
respectively and results are shown for a double powerlaw LF with various values of a and AM. 



the 1, 2, 3 and 4 a levels of a Gaussian distribution. The up- 
per two and lower two rows of figure 5 correspond to velocity 
and density evolution respectively and in each case results 
are shown for a double powerlaw LF with various values 
of q and AM. The figure shows that the preferred values 
are found near near 7 = and (3 = —2.1. We also show 
contours of the multiple image fraction (light contours); the 
solid, dashed, dot-dashed and dotted contours correspond 
to lens fractions of 0.1, 0.03, 0.01 and 0.003. We find that 
the multiple image fraction should be ~ 1 — 3%. This value 
is lower than previous estimates due to constraints on the 
shallow luminosity function, and unfortunately implies that 



az~6 lens may not be found among the complete sample 
of SDSS z ~ 6 quasars. 



8 WHAT IF SDSS J 1148-5251 WERE 
MULTIPLY IMAGED? 

While current high resolution imaging suggests that all of 
the z ~ 6 quasars are point sources, White et al. (2003) 
note the possibility that SDSS J1148-5251 is multiply im- 
aged cannot be ruled out by current observations due to the 
small angular diameter of the Einstein ring in this system; 
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Figure 6. Differential probability distributions for f3. The lensing constraints that include information on the alignments between singly 
imaged quasars and foreground galaxies as well as a single multiply imaged quasar are denoted by the dark lines. The solid grey curves 
correspond to the probability functions for /3 based on the luminosity distribution alone (Fan et al. 2003). Results arc shown for a double 
powerlaw LF for various values of a and AM. The upper and lower two rows show results assuming density evolution with 7 = and 
7 = 2 respectively. 



A8 ~ 0.3" for a a = 250km/sec galaxy at z = 4.94 lensing a 
quasar at z = 6.37. For comparison the z ~ 6 quasars have 
been imaged at a resolution of ~ 0.4" (Fan et al. 2003). It 
is therefore possible that SDSS J1148-5251 is multiply im- 
aged but appears as a point source, though we note that 
this is an unlikely scenario since with 4" seeing, a double 
with a 3" separation should be recognisable if the flux ratio 
is smaller than 10 : 1 (Chris Kochanek 2003, private com- 
munication). Multiple imaging of SDSS J1148-5251 would 
have important implications for the study of the z ~ 6 LF, 
invalidating the constraints on /3 obtained in § 5. We have 



therefore computed the limits imposed on /3 by the observa- 
tion of one multiply imaged source (SDSS J1148-5251) and 
one singly imaged source with high alignment (SDSS J1044- 
0125) among a sample of six quasars at z ~ 6. 

The likelihood function for this scenario may be writ- 
ten as the product of the probability that a source is lensed 
with the likelihood that the remainder are singly imaged. 
The observation that one of the quasars is multiply imaged 
constrains large values of /3, while the observation of high 
alignment without multiple imaging limits small values of /3 
as discussed in previous sections. The maximum of the com- 
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bined likelihood function lies in the lower end of the range 
specified by the luminosity distribution (Fan et al. 2003). 
The posterior probability function for the sample of z ~ 6 
quasars may be written 



dP 

1(3 = Nexp 




xpi cns (/3, Ll) X 
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(27) 
(28) 



and is plotted in figure 6 for various values of AM and a. We 
find the preferred value in this case to be (3 ~ —3. As noted 
in § 2.1 it is surprising that one quasar in the sample would 
be lensed, but not at a highly significant level. The obser- 
vation of one lensed lensed quasar therefore prefers smaller 
values of (3, for which the magnification bias is larger, and 
also slightly tightens the allowed range for j3. In summary, if 
SDSS J1148-5251 were multiply imaged the preferred value 
for the slope would be (3 ~ —3, which is ruled out at the 
90% level if the quasar is singly imaged but with a high 
alignment. From figure 5 we see that (3 ~ — 3 implies a lens 
fraction of ~ 3 — 10% rather than ~ 0.3 — 1%, which is 
more consistent with previous estimates. This underlines the 
importance of determining whether this quasar is multiply 
imaged or merely lensed. 



9 CONCLUSION 

From their analysis of the luminosity distribution of quasars 
at z ~ 6 Fan et al. (2003) determined a slope for the quasar 
LF of (3 i; —3.3 (90%). This slope is consistent with the value 
found for the slope of the LF at z ~ 4.3, but is not consistent 
with the slope of the LF of bright quasars at z 5s 3. It is also 
possible to constrain the slope of the LF using the fraction 
of multiply imaged lensed quasars. Fan et al. (2003) com- 
puted the probability of obtaining a lens fraction of zero as 
a function of f3 and found that the constraints were similar 
to those of the LF. We have performed a Bayesian analy- 
sis including both of these (nearly) independent constraints, 
yielding the result that at 90% confidence (3 J> —3.3 — > —3.0 
provided that the quasars are at least 2 magnitudes brighter 
than the unknown position of a break in a double powerlaw 
LF. The systematic dependence in the constraint is due to 
the unknown slope of the LF at fainter luminosities, the lu- 
minosity of the break and the uncertain evolution in the lens 
galaxy population. 

While inclusion of constraints from the multiple image 
fraction somewhat improves the limits on (3, we have shown 
that the additional information from observations that nei- 
ther of the two quasars which lie near to the line of sight 
to foreground galaxies (and which are therefore lensed) are 
multiply imaged provides a stronger lensing based constraint 
on the slope (3. We find that for a double powerlaw LF 
(3 > -3.1 -» -2.9 with 90% confidence. Unlike the con- 
straint that uses only the multiple image fraction, this limit 
is nearly independent of evolution in the lens population, 
and adds further evidence of a trend to shallower LF slopes 
at large redshifts. 



We also find that the existence of a lens galaxy at z ~ 5 
in a sample of two lenses constrains the evolution in the 
massive galaxy population to be close to that of constant co- 
moving density (provided that the quasar behind the z ~ 5 
galaxy is subject to magnification bias). This lack of evo- 
lution in the lens population is consistent with the lack of 
redshift evolution in the velocity function of dark-matter ha- 
los (for velocity dispersions near 200km sec -1 ) as predicted 
by the Press- Schechter formalism. 

Finally the constraints on the quasar luminosity func- 
tion and lens population have been used to compute an im- 
proved estimate for the expected z ~ 6 multiple image lens 
fraction of ~ 1 — 3%. This value is lower than previous esti- 
mates due to the tight constraints on the slope of the LF. 
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